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Abstract 

We provide a rather complete description of the sharp regularity theory for a family of 
heterogeneous, two-phase variational free boundary problems, — > min, ruled by nonlin- 
ear, degenerate elliptic operators. Included in such family are heterogeneous jets and cavities 
problems of Prandtl-Batchelor type, 7=0; singular degenerate elliptic equations, < 7 < 1; 
and obstacle type systems, 7=1. Linear versions of these problems have been subjects of 
intense research for the past four decades or so. The nonlinear counterparts treated in this 
present work introduce substantial new difficulties since most of the classical theories devel- 
oped earlier, such as monotonicity and almost monotonicity formulae, are no longer available. 
Nonetheless, the innovative solutions designed in this article provide new answers even in the 
classical context of linear, nondegenerate equations. 
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1 Introduction 

Let a C M" be a bounded domain, 2 < p < +00, f e L«(^2) forq>n and (p G nL°°(f2), 
with, say, (p^ ^ 0. The objective of the present manuscript is to derive optimal interior regularity 
estimates for the archetypal class of heterogeneous non-differentiable functionals 

(1.1) A(^) := 1^ (|Vvr + Fy(v) +/(X) ■ v) dX min, 

among competing functions v G Wq'^{Q.) + 9. The parameter 7 in (11.11) varies continuously from 
to 1, i.e., 7 G [0, 1] and the non-differentiable potential Fy is given by 

(1.2) Fy{v):=X+{v+y + X-{v-y, 

for scalars < X- < X+ < °°. As usual, := max{±v,0}, and, by convention, 

(1.3) Fo{v) := KX{v>o}+^-X{v<o}- 

The non-differentiability of the potential Fy impels the Euler-Lagrange equation associated to 
to be singular along the a priori unknown interface 

57:= (^{My> 0}U(9{My < 0}) na, 

between the positive and negative phases of a minimum. In fact, a minimizer satisfies, in some 
weak sense, the following p-degenerate and singular PDE 

where ApU denotes the classical j!?-Laplacian operator, 

Apu:=di\{\Vu\P-^Vu). 
The potential Fq is actually discontinuous and that further enforces the flux balance 

(1.5) \Vu+\P - I 1'^ = ^ (A+ - A_) , 

p — l 

along the free boundary of the problem, which breaks down the continuity of the gradient through 
do. 

A number of important mathematical physics problems, coming from several different con- 
texts, are modeled by optimization setups, for which equation (11.11) serves as an emblematic, lead- 
ing prototype. This fact has fostered massive investigations, and linear versions, p = 2, of the 
minimization problem (11.11) have indeed received overwhelming attention in the past four decades. 
The upper case 7 = 1 is related to obstacle type problems. The linear, homogeneous, one phase 
obstacle problem, i.e., p = 2, f{X) = and (p >0 was fully studied in the 70's by a number of 
leading mathematicians: Frehse, Stampacchia, Kinderlehrer, Brezis, Caffarelli, among others. It 
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has been established that the minimum is locally of class C ' and this is the optimal regularity for 
solution. The two-phase version of the problem, i.e., with no sign constrain on the boundary datum 
(p, challenged the community for over three decades. C^'^ estimate for two-phase obstacle prob- 
lems was established in GOll with the aid of the powerful almost monotonicity formula obtained in 

m. 

The lower limiting case, 7=0, relates to jets flow and cavities problems. The linear, homo- 
geneous, one phase version of the problem was studied in [[T]|, where it is proven that minima are 
Lipschitz continuous. The two-phase version of this problem brings major new difficulties and 
C'^-^ local regularity of minima was proven in [2J, with the aid of the revolutionary Alt-Caffarelli- 
Friedman monotonicity formula, developed in that very same article. Gradient estimates for two- 
phase cavitation type problem with bounded non-homogeneity, i.e., p = 2, f E U°, 7 = in (II. 1|) . 
was established by Caffarelli, Jerison and Kenig with the aid of their powerful almost monotonicity 
formula, jSl. 

The intermediary problem < 7 < 1 has also received great attention in the past decades. The 
related free boundary problem can be used, for example, to model the density of certain chemical 
specie, in reaction with a porous catalyst pellet. The linear, p = 2, one-phase, <p > 0, homogeneous, 
/ = 0, version of the problem (|l.ll) is the theme of a successful program developed in the 80's by 
Phillips and Alt-Phillips, lfT9l . ifTSl and |l3l, among others. In similar setting. Holder continuity 
of the gradient of minimizers was proven Giaquinta and Giusti ifTOl . Further investigations on 
the linear, two-phase version of this problem also require powerful monotonicity formulae in their 
studies, see ETll . 

In the mathematical analysis of variational free boundary problems as (11.11) . the first major key 
issue to be addressed concerns the optimal regularity estimate available for a given minimum. A 
simple inference on the weak Euler-Lagrange equation satisfied by a minimum. Equation (11.41) and 
also the Flux Balance (11.51) for 7=0, revel that ApU blows-up along the free boundary of the prob- 
lem, •= d{uY > 0} U d{uY < 0}. Therefore, it becomes a fundamental question to understand 
precisely how this phenomenon affects the (lack of) smoothness properties of minima. Under such 
perspective, and to some extent, the theory of two-phase free boundary problems governed by 
non-linear, degenerate elliptic operators had hitherto been unaccessible through current literature, 
mainly due to the lack of monotonicity formulae in this context. 

In the study of sharp smoothness properties of minima to the functional further difficulties 
also arise from the very complexity of the regularity theory for the governing operator Ap. We 
recall that p-harmonic functions, i.e., solutions to the homogeneous equation 

Ap/z = in^i, 

are locally of class for an exponent < Up < I that depends only upon dimension and 

p. The precise value of Up is in general unknown - see [[T4l| for the planar case n = 2. This 
fact indicates that interior estimates available for p-harmonic functions, that in turn are below 
quadratic, ^ will compete with optimal growth along the free interface 5^y. The regularity theory 
for heterogeneous equations Ap^ = f{X) is even further involved and, up to our knowledge, the 
understanding on this class of problems is not yet fully complete. 
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The fundamental contribution of the present work relies on a complete description of the under- 
lying sharp regularity theory for minima of which primely opens the possibility to investigate 
these classes of problems in full generality. 

From the mathematical point of view, the exponent 7 appearing in (11.11) should be compre- 
hended as the parameter that measures the singularity of the absorption term of the related equa- 
tion. For non-differentiable but continuous functionals, with < 7 < 1 , it has been conjectured 
that the gradient of a minimum is locally Holder continuous, even through the singular free inter- 
face g^y. The first result we present in this paper gives an affirmative answer to such question. 
Furthermore, it provides the asymptotically optimal C^'" interior regularity theory available for 
minima of such functionals. 

Theorem 1.1 (C^'" regularity estimates). Let u he a minimizer of the prohlem (11.11) . Assume 
< 7 < 1 and f G L^(f2), for some q > n. Then u G c]^'^, for 

7 {q-n) 



(1.6) a := < ttp 



Furthermore, for any Ol d O., there exists a constant C > depending only on, Q.', n, p, q, 
ll9llL~(a). ||/||L9(a). and {Up - a), such that 



<C. 



Let us further explain that the asymptotically sharp regularity estimate indicated in (11.61) should 
be read as 



(1.7) 



If min^^,^}<a,„ then ^GC^HA'aHi}. 
If min I I > Up, then u G C^'^, for any < a < Up. 



Before continuing, let us make few comments on Theorem 1 1.1 1 and its implications. The key 
ingredient of the regularity estimate established in Theorem 11.11 reveals how the competing forces 
involved in the lack of smoothness for minima of (11.11) . namely 

(regularity theory for A^) x (singular absorption term ~ w^^^) x (roughness of the source /) 

get adjusted, via the sharp relation (11.61) . Regarding the exponent Up, one easily verifies that for 
some constant Op, the function X dp{X^)~p^ is p-harmonic, thus oCp is never beyond the critical 
value ^ri. However, nonnegative functions with bounded p-Laplacian, v, do grow as dist'^ i (X, 
away from = d{v > Q}, see [[TT|. In this particular setting, it is possible to replace ap by 

y 

in (|1.6I) . Thus, at least if / G L°°, Theorem 11.11 revels u G Cp-^ , which is the precise generalization 
of the optimal regularity estimate obtained for the one-phase linear setting p = 2, see for instance 

imiiiii. 
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Confronting the effect of the singular absorption term ~ ^ and the influence of integrabil- 
ity properties of the source /, we conclude that solutions to (11.11) are locally in (7^'^"^"^ 
provided f E for any 

(1.8) g>n.^^=:q{p,n,r). 
Interestingly enough, one verifies that 

{q{p,n,r)-n)p 1 



(1.9) q{p,n,\ ) = oo and lim 



y^\{p-l)q{p,n,y) p-V 
Also it is revealing to compute the limit 
(1.10) \m\q{p,n,'Y) = n, 

which leads us to the discussion of the delicate limiting case, 7 = in the minimization prob- 
lem (|l.ll) . As mentioned earlier in this Introduction, for homogeneous, / = 0, linear, p = 2, jets 
and cavities problems, Lipschitz regularity estimates have been established in the one-phase and 
two-phase case, respectively in lUl and [|2||. Heterogeneous, two-phase versions of the problem 
could only be approached quite recently, with the aid of the almost monotonicity formula, [[5]|. 
However, the Caffarelli-Jerison-Kenig monotonicity formula requires a one-side bound for the 
non-homogeneous term f{X), namely, f{X) > — C. Thus, even for linear problems, p = 2, Lip- 
schitz estimates for minimizers of (11.11) . 7=0, are only known if / G L°°(t2). We further point 
out that the integrability exponent obtained in (11.101) is a borderline condition, as it divides the 
regularity theory for (non-singular) Poisson equations, Lu = f, between continuity estimates when 
/ G L"^^ and differentiability properties when / G L"+^. The optimal regularity theory for the con- 
formal case / G L" is rather delicate. It has been recently established by the third author, [|22l . that 
solutions to nonlinear equations F{X,D^u) = f{X) G L" has a universal Log-Lipschitz modulus of 
continuity, i.e., 

\uiX)-uiY)\<\X-Y\-\og\X-Y\. 

Such regularity is optimal in the context of heterogeneous equations with L" right-hand-sides. 
After some heuristic inferences, it becomes reasonable to inquire whether minimizers of problem 
(11.11) . with 7=0, also has a universal Log-Lipschitz modulus of continuity. The second main result 
we establish in this paper states that indeed minimizers of ^0 with sources / G L" also enjoy such 
an optimal universal modulus of continuity. 

Theorem 1.2 (Log-Lipschitz regularity for 7 = 0). Let u be a minimizer of the problem (11.11) . with 
7=0 and f G L"{Q). Then u is Log-Lipschitz continuous and for any Q! d Q., there exists a 
constant C that depends only on, Q.', n, p, ll/llL"(n)' ^+ ^^'^^ ^^'^^ 

\u{x)-u{Y)\ <c|x-y|iog|x-y|. 

In particular. Theorem 11.21 assures that u G c{^^(f2) for any T < 1. We further mention that 
Theorem 11.21 is sharp due to the borderline integrability condition on the source /. It becomes. 
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therefore, a major issue to investigate whether functional has a locally Lipschitz minimizer, 
provided / G L'^{Q.) for q > n. We highlight that this question remained open even for the linear 
case p = 2, as no almost monotonicity formula can be established unless the source function / is 
bounded. 

The final ultimate goal of the present work gives an affirmative answer to the above question. 

Theorem 1.3 (Existence of Lipschitz minimizer for 7 = 0). Assume f E LF{^), q> n. Then, there 
exists a minimum uq to the functional that is locally Lipschitz continuous. Furthermore, for 
any Q! (e Q., there exists a constant C > 0, depending only on, O! , n, p, |<p||L~(f2). 
and X-, such that 



We leave open the question whether any minima of the functional is Lipschitz continuous. 
A critical analysis on the machinery employed in the proof of Log-Lipschitz estimates, Theorem 
11.21 reveals that it should not be possible to access the CP'^ regularity theory for minima of 
through pure energy considerations, even if the source f E LT. The decisive, innovative path sug- 
gested in this work to prove Theorem 11.31 uses a new varying singularity perturbation technique, 
cp^r ~^ c^Oj which, in turn, relies on geometric non-degeneracy properties of the limiting mini- 
mizer. 

To finish up the Introduction we mention that the Theorems established in this article extend, 
with minor modifications, to further involved energy functionals of the type 



where G is a p-degenerate kernel with C coefficients, \Gy\ < Fy and v)| < g{X)\v\'^, where 
<m < p and g E L^, for q > max | ^^z^^^^j- We have chosen to present our results in a simpler 
setting as to further emphasize the novelties and core ideas of the methods and techniques designed 
in this work. 

The paper is organized as follows. In Section [2] we gather few tools that we shall use in the 
proofs of Theorem 11.11 Theorem 11.21 and Theorem II. 3[ In Section [3] we comment on existence 
and establish universal L°° bounds for minima of problem (11.11) . Section |4] is devoted to the proof 
of Theorem 1 1.1 1 and in Section [5] we establish Log-Lip estimates for cavitation problems, proving 
therefore Theorem 11.21 Under the condition / E L^, q > n, in Section [6] we show sharp linear 
growth and strong nondegeneracy properties for solutions to the cavitation problem 7 = 0. In Sec- 
tion|7]we investigate stability properties for the family of free problems in terms of the singular 
parameter < 7 < 1 . More precisely we show that local minima of functional converges to 
a local minima of the functional ^0, as 7 — )■ 0. In the final Section [8] we conclude the proof of 
Theorem 1 1.3 1 

Acknowledgement RL research has been funded by Capes and CNPq. OSQ has been partially 
supported by FAPESP/SP-Brazil. ET thanks support from CNPq-Brazil. 
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2 Preliminaries and some known tools 

In this section we gather some preliminaries results that we will systematically use along the arti- 
cle. Initially, as mentioned within the Introduction, clearly one should not expect solutions to the 
minimization problem (11.11 ) to be smoother than p-harmonic functions. Therefore, the regularity 
theory for degenerate elliptic operators is a first key ingredient in understanding sharp estimates 
for minima of ^j. 

There are several different strategies to establish the C^'"/' regularity theory for p-harmonic 
functions, see for instance Q, |[8l, [fTSl , l|23l , l|25l and [|26ll . We state such result for future refer- 
ences. 

Theorem 2.1 (C^'" estimates for p-harmonic functions). Let h G W^'P{Bi) satisfy Aph = in Bi 
in the distributional sense. Then, there exist constants C > and < < 1, both depending only 
on dimension and p, such that 

ll^llcl'«/^(Si/2) - ^ll^llz^(si)- 

A particularly interesting approach was suggested by Lieberman in lfT6l . where the regulari- 
ty theory for p-harmonic functions is accessed through the following leading integral oscillation 
decay lemma: 

Lemma 2.2 (Lieberman, [|16,1 . Lemma 5.1). Let hbe a p-harmonic function in Br C M". Then, for 
some positive constant < (Xp < I, there holds 

r / r\ n+pa,, r 

\Vh{X) - {Vh)/dX < C (-) \Vh{X) - {Vh)^\PdX, 

where C = C{n,p) > is a positive constant. 

In Lemma [Z2] and throughout this article we use the classical average notation 

Local Holder continuity for heterogeneous equations Ap^ = f can be delivered by means of 
Hamack inequality, which will be another fundamental tool in our analysis. 

Theorem 2.3 (Hamack inequality, e.g. (HII). Let ^ G W^-p{Br), ^ >0 a.e., satisfy ApE, = f in 
Br in the distributional sense, with f & L^ (Br) and Q > jy Then, there exists a constant C, > 
depending only on n, q, p and R — r such that 




for allO<r<R. 

In the sequel, let us discuss some further inequalities that will be used in the proofs of our 
main results. The estimates presented herein have elementary character and are mostly known. We 
include them for completeness purposes and courtesy to the readers. 
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Lemma 2.4. Let \]/ e with 2< p and h e W^'^{Bi) solution to Aph = in Bi. Then, for 

a constant c = c{n,p) > 0, there holds 

I [\V\if\P-\Vh\P]dX>c [ \V{\ir-h)\PdX. 

Proof. For each < T < 1, let 0^ denote the linear interpolation between i// and h, i.e., 1//^ : = 
(1 — x)h. From Fundamental Theorem of Calculus we have 

(2.1) / {\V\if\P-\Vh\P)dX= —( f \V\ifr\PdX^dT. 

Jbi Jo dx \Jbi J 

Passing the derivative through and using the fact that div (| V/z|p^^V/z) ■{\ir — h) = in 5i, we find 
(2.2) 



lllxi^j I^V^^I''^^)^^ = pJq^^J {\^¥^\'''^^¥T-\^h\P-^Vh)-y{\ir-h)dX 

= pC-dxl {\V\ir,\P-^V\ir,-\Vh\P-^Vh)-V{\if,-h)dX, 

Jo T Jb, 



because — h = x{^f — h) . The Lemma now follows easily from the well known classical mono- 
tonicity 

(2.3) m'-^^x - \^2r%. ^1 - ^i) > c{n,p)\^i - ^2^, 

for any pair of vectors i^i, (^2 £ I^"- In fact, combining (12.11) . (12.21) and (12.31) we reach 

/ {V\if\P-\yh\P)>c{n,p) [\-Ut [ \V{\ir^-h)\P = c{n,p) f\p-UT [ |V(v/-/z)|^, 
Jbi Jo Jbi Jo Jbi 

and the Lemma follows. □ 
Lemma 2.5. Let 7 G (0, 1). For any positive scalars a> 0, b > there holds 

(2.4) {a + by <{a^ + b^). 
Proof. In fact, just notice that, since 7— 1 is negative, we have 

(2.5) t^-^>{t + ay-\ V?G(0,oo). 



Then, integrate (12.51) from to to obtain the desired inequality. □ 

Next we prove two useful asymptotic inequalities. 
Lemma 2.6. Let < /i < 1 and suppose a real function verifies 

<^)(r) <A(/>0(r)^ + r"2), 



for r small enough. Then (r) = O ( j-™"!,'^^' 1 1 ^ approaches zero. 
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Proof. In fact, if ^(r) < r^i0(r)^ + then for j8 := min {e2, : there holds 



< 



< 



+ 1 



+ 1, 



since < /3, The above readily implies 0(r) = O(r^) as claimed. 

Lemma 2.7. Le? ^{s) be a non-negative and non-decreasing function. Suppose that 



□ 



(2.6) 



(^)(r)<Ci 



LVi? 



for all r < R < Rq, with Ci, a,/3 positive constants and C2,/i non-negative constants. Then, for 
any o < min{a,/3}, there exists a constant jUq = {Ci, a,fi,o) such that if ji < IM), then for all 
r < R < Rq we have 



(2.7) (r) < C3 j [<^)(7?)+C2i?^] 

where C3 = C3 (Ci , o — min{a, /3 }) a positive constant. In turn, 

(2.8) <^)W<C4r^, 
where C4 = €^{€2, C^^Rq^^^o) is a positive constant. 

Proof. We can assume /3 < a and, in this case, it suffices to show the estimate for C7 = /3. For 
< < 1 and R<Rov/e have 



(2.9) 



^{GR) < Ci 



7? 



^{R)+C2R^ 



0"Ci [1+Ai0""] 0(i?)+C2i?^. 



We choose < < 1 such that 2Ci0" = 0^ with j5 < 5 < a. Now we take /io > satisfying 
jUo0~" < 1- Thus we obtain for all ^ < ^0 

(2.10) (j){eR) < d^^{R)+C2R^. 
Inductively we get 

(2.11) ^{d'^^'^R^ < d^^{d^R^+C2d^^R^ 



< C3d^''+^^^\(j){R)+BRP 
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for all keN. Hence, taking k such that O'^+^R <r< O'^R we obtain (lO) . 
Finally, we have 




which proves inequality (|2.8I) . 



□ 



3 Existence and L°° bounds of minimizers 

In this section we establish existence and pointwise bounds for a minimum of the functional ^y. 
The arguments presented herein works indistinctly for the cases < 7 < 1 and 7=0. 

Theorem 3.1 (Existence and LT bounds). Let Q. C M" be a bounded domain, f E L^{^), q > n, 
cp E W^'P (Q,) nL°° (Q.) and < 7^ A_ < 00 be fixed. For each < 7 < 1, there exists a minimizer 
Uy to the energy fiinctional 



Initially we show that Iq > —0°. Indeed, for any v eWq'^ + (p, by Poincare inequality there exists 
a positive constant c = c{n,p,Q., \\f\\ui) > such that 





Proof. Let us label 




(3.1) 



c V 



p 

LP 



c||<l>||£,.-||V0||f,<||Vv||£ 



By Holder inequality, since 



(3.2) 



q>n> 



P 



we have 



/ f{X)vdX < 11/11 



LP 



p \\v\\lp < Ci{n,p,Q.)\\f\\Lq\\v\\LP. 



REGULARITY FOR TWO-PHASE FREE BOUNDARY PROBLEMS 11 

Young's inequality yields the existence of a positive constant C = C{n,p, Q., WfWui) > such that 

(3.3) -c\\v\\l,,-C <C,{n,p,a)\\f\\t4v\\iJ'- 
Combining (13.11) and (13.31) we obtain 

(3.4) _C-c||<|)||£,-||V<|)||f,<||Vv||^.-Ci(n,p,f2)||/|U,||v||z;^ 
which reveals 

-C-c||<|)||f,,-||V(^)||^, < ||Vv||^,-Ci(n,;,,^2)||/||^.l|v||LP. 

Finally, we reach 

(3.5) -C-cWni^ - \\W(^\\Pj, < \\Vv\\Pj,-C,{n,p,mf\\L4ALr < /M- 

Let us now show existence of a minimum. Let Vj G W^'^{Q.) be a minimizing sequence. For 
7> 1, 

M^i) < ^0 + 1. 

From (13.51) and Holder inequality we obtain 



(3.6) 



l^\Vvj\PdX < C\\vj\\u^+Io + l. 



By Poincare inequality we estimate 

(3.7) C\\vj\\lp < C(||Vv,||l. + ||V(|)||,;. + ||(^||z;.). 
Also we have 

(3.8) C||VvJz,p <C + ^||Vv^-||£,. 
Combining (Il6l) . (13771) and (|3^ we reach 

(3.9) / \Vvj\PdX < C{\\V(^\\u, + U\\lp)+Iq+1. 

Ja 

Thus, using Poincare inequality once more, we conclude that {vj — ^} is a bounded sequence in 
Wq'P (Q.). By reflexivity, there is a function u G W^'^ (Q.) such that, up to a subsequence, 

Vj^u weakly in VK^'^ (^2) , vj^u inL''(f2), vj^u a.e.inQ,. 
From lower semicontinuity of norms, we readily obtain 

/ \Vu\PdX < liminf / |Vv,|W. 
Ja Jo. 
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By pointwise convergence we have, in the case < 7 < 1 , 

/ Fy(u) +f(X)udX < liminf I Fyivi) +f(X)vjdX. 
Jq. Ja 

For 7=0, recalling that we are working under the regime A+ > A_, we have. 



Thus, 



J^^-X{.<o}dX < liminf (/ 



\J{u<0} 

On the other hand, since Vj^u a. e. in we have 



lim / A+^r 



Hence, 



In conclusion. 



/ Fo{u)dX < liminf / Fo{vj)dX. 
JQ. Jo. 



/y{u) < liminf ^y(v^) =/o, 

for < 7 < 1 , which proves the existence a minimizer. 

Let us now turn our attention to L°° bounds of Uy, which hereafter in this proof we will only 
refer as u. Let us label 



7o ■-- 



sup0 

do. 



that is, the smallest natural number above sup^^^ 0. For each j > jo we define the truncated function 



(3.10) Uj 



j-sing(M) if |m|>7 

u a \u\ < j. 



where sing(M) = 1 if w > and sing(M) = — 1 else. If we denote := {\u\ > j}, we have, for each 
j > jo 

(3.11) u = Uj in and Uj = j-sing(M) in Aj. 
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Thus, by minimality of u, there holds, for < 7 < 1, 
(3.12) / \Wu\PdX = I \Wu\P -\Wui\PdX 

JAj Jo. 

< I f{uj-u)dx+ I x+ ({u+y - {u+)A dx 

JAj JAj 

+ X-(^{ujy-{u-yyx. 



Notice that 



/ f{uj-u)dX = [ f{j-u)dX+[ f{u-j)dX 

JAj JAjn{ii>o} JAjn{u<0} 

< if \f\i\u\-j)dX. 

JAj 



Moreover, we have 



X+ [ ({uiy-{u+)AdX = A+ / {f-\u\^)dX 

JAj V J J >/Ajn{i(>0} 

+ A+/ [i^-jY)y-{u^)y)dx 

JAj^{u<X)) 



< 



and 



xf ({ujy-{u-y)dx = xf {uyy-{u-y)dx 

JAj V J J >/Ayn{l(>0} 

+ A- / {p-\u\^)dX 
JAjr\{u<sS\ 

< 0. 

Then, we find 

(3.13) / Fy{Uj)-Fy{u)<0. 

JAj 

For 7 = it suffices to notice that Uj> and u have the same sign. From the range of truncation we 
consider, it follows that — j)^ G Wq'^ (Q). Hence, applying Holder inequality and Gagliardo- 
Nirenberg inequality, we find 



[ \f\{\u\-jydx < 11/11 .iKi^i-j) 

Jaj 
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where p* := Young inequality gives, 



1 



(3.14) ||/||l.|A/ -* 5||Vi.||^,^^.)<C|A,|'-i-^(^-i)-W.-i)+-||Vi.||^,^^^^. 
Combining (13.121) and (13.141) we obtain 

(3.15) / \Vu\PdX <C\Aj\^-"+', 

Jaj 

where e = and (see (lO) and (Il9l) substituting /q by ^7(^>)) 

(3-16) """i'('4,,o) - I^./o|'^II"IIl/'(a,.J 

Boundedness of m now follows from a general machinery, see for instance, ll24l . Chap. 2, Lemma 
5.2, Page 71. □ 

Remark 3.2. A consequence of L°° estimates for a minimum u to the functional the universal 
control of u in W^^P. In fact, we have 

(3.17) [ \Vu\Pdx < Jy{(p)- f Fy{u)dX+ I \f{X)\\u\dX 
Jo. Ja Ja 

< M'P)+C{n,PM\f\\L'>) 

< c, 

where C = C(n, p, (p, ||/||/^9) > is a positive constant. Here we used the elementary inequality 
< max { 1 , ? }, for ? > and < 7 < 1 . In conclusion, 

(3.18) \\u\\wi.p<C. 

We close up this Section by stating the Euler-Lagrange equation associated to the functional 
^y, < 7 < 1 as well as the flux balance - also known as the free boundary condition - satisfied 
by a minimum uq to ^0, through the free boundary. The proofs of this facts are rather standard 
and we omit them here. 

Proposition 3.3. Let Uj be a minimum to the functional ^y, < 7 < 1. Then Uy solves 

(3.19) ^^u = ^{x+{u^y-'x{u>o)-^-{u-T'x{um) + \f^^) 

in the distributional sense. Also, ifuo is a minimum of ^q, with \ {uo = 0}| = 0, / G (t2), q> n, 
Xq G ^^{uq) U ^~uq) a generic free boundary point and B a ball centered at Xq. Then for any 
$GC(J(5,R"), there holds 

lim / (((p-l)|Vwor-A+)vi,$)J^"-i 

£l\0 J 

Sn{Ko=ei} 

+ 

lim / {{{p-l)\VuQ\P-X.)v2,^)dje''-' 
Bn{«o=e2} 

= 0, 
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where Vi and V2 denote the outward normal vector on 5 fl {mq = £1 } (^nd 5 fl {mq = £2} respectively. 
In particular, the flux balance 

\Vu+\P-\Vu,\P = ^iX+-X.), 
holds along any C^'" piece of the free boundary. 

4 Sharp C^^ estimates for minima 

This Section is devoted to the proof of Theorem 11.11 which assures optimal Holder continuity 
estimates for the gradient of minima of the energy functional ^y, for < 7 < 1 and q > n. The 
borderline situation 7 = and / G L" will be addressed in the next Section. 

Hereafter in this Section, u = Uy denotes a minimizer of the functional ^y, with < 7 < 1. 
Theorem I 1 . 1 [ concerns an optimal interior regularity result; therefore, in order to prove such interior 
estimate, we fix an arbitrary point G t2 and R>Q such that R < dist(Xo, dQ.) . We will show that 
u G C^^" at Xq, for a as in (fTel) . 

In the sequel we show the first main step in our strategy to obtain sharp regularity theory for 
minima of the energy ^y. 

Lemma 4.1 (Comparison with p-harmonic functions). Let u G W^'P{Bji) and h G W''^(5/e) satisfy 
Aph = in Br in the distributional sense. Then, there exists a positive constant C = C{n,p) > 
depending on dimension and p such that for each < r <R, there holds 

(4.1) / \Vu{X)-{Vu)\PdX < C(-) / \Vu{X)-{Vu)JPdX 

J Br ^R'' J Br 

+ C [ \Vu{X)-Vh{X)\PdX, 

JBr 

where Q < dp < I is the optimal exponent in Lemma \2?2\ which, in turn, reveals the sharp C^'" 
estimate stated in Theorem \2.1\ 

Proof. For each r G (0,i?] we estimate, 

(4.2) / \Vu{X)-{Vu)fdX < Cp [ \Vu{X)-{Vh)fdX 

JBr JBr 

+ Cp[ \iyu)^-{vh)fdx, 

JBr 

for a constant Cp that depends only on p. Analogously, we obtain 

(4.3) f \Vu{X)-{Vh)ydX < Cp f \Vu{X)-Vh{X)\PdX 

J Sf J Bf 

+ Cp f \yh{x)-{yh)fdx. 

JBr 
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In the sequel, we apply Holder inequality and estimate 

1 



(4.4) / \{Wu),.-{Wh)XdX 



\B, 



< 



< 



1 



Br 



{Vu{X)-Vh{X))dX 



Br 



\yu{X)-Vh{X)\dX 



\BAP-' ] \JBr 



\Vu{X)-Vh{X)fdX 



L-\ P 

p 



Br 



\Vu{X)-Vh{X)fdX. 



Combining (|4!2l) . (|431) and (l4i4l) we obtain 

(4.5) [ \Vu{X)-{Vu)fdX < Cp I \Wh{X)-{Wh)fdX 

+ Cp [ \Vu{X)-Vh{X)\PdX. 

JBr 

Interplaying the roles of u and h in (14.51) and arguing in the bigger ball Br, we find 

(4.6) f \Vh{X)-{Vh)^\PdX < Cp j \Wu{X)-{Wu)^\PdX 

JBr JBr 

+ Cp [ \Vu{X)-Vh{X)\PdX. 

JBr 



Now, in view of Lemma 12 . 21 and (14.51) we can further estimate 

r \ n+pap 



(4.7) 



/ \Vu{X)-{Vu)\PdX < C{n,p)[-] / \Vh{X)-{Vh)r,\PdX 

JBr Vi?/ J Br 

+ C{n,p) [ \Vu{X)-Vh{X)\PdX. 

JBr 



Hence, combining (14.61 ) and ( 14.71 ) we readily obtain 

(4.8) / \Vu{X)-(yu)ydX < C{n,p){-) / \Vu{X) - {Vu)JP dX 

JB. V^/ JBr 



which finally implies 



+ C{n,p) 



Br 



Br 



\Vu{X)-Vh{X)\PdX, 



r / r\ n+pa,, r 

(4.9) / \Vu{X)-{Vu)fdX < C(-) / \Vu{X)-{Vu)j^\PdX 

JBr V/v/ J Br 



+ C [ \Vu{X)-Vh{X)\PdX, 

JBr 



and the proof of Lemma l4.1l is concluded. 
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We have now gathered all the tools and ingredients we need to establish local Holder continuity 
of the gradient of a minimum of the energy functional < 7 < 1 . 

Proof of Theorem We start off the proof by denoting, for writing convenience, Br := Br{Xq) 
and M = My a given minimum of the functional < 7 < 1 . Let h be the p-harmonic function 
in Br that agrees with u on the boundary, i.e., 

Aph = OmBR and /? - m e '^(5^). 

By Lemma l4n we have 

r / r\ n+pa,, r 

(4.10) / \Vu{X)-{Vu)ydX < C -) / \Vu{X)-{Vu)j.\PdX 

J Br ^R'' J Br 

+ C [ \Vu{X)-Vh{X)\PdX. 

J Br 

On the other hand, by the minimality of u we have 
(4.11) 



f {\Vu\P-\Vh\P)dX< f {Fy{h)-Fy{u))dX+ f f{X){h-u)dX. 

J Br J Br J Br 

Invoking Lemma [Z4l there exists a constant C3 = C2,{p,n) > such that 
(4.12) C3 / (\Vu\P-\Vh\P)dX> [ \V{u-h)\PdX. 

J Br J Br 



Moreover, we have 



Br 



Fy{h) - Fy{u)dX = X+ I {h+) ^ - ^ 



Br 



dX + X- 



Br 



{h-y-{u-y 



dX 



with 
(4.13) 



Br L 



dX < 



+ 



{h+>u+} 
{h+>u+}n{u+>0} 
{/i+>t(+}n{;(+=0} 



dX 



{h+-u+ydX. 



dX 



Notice furthermore that 
(4.14) 



'{h+>u+}n{u+=o} 
By Lemma [231 there holds 



(h+-u+ydX< [ ih-u) 

J{h+>u+}r\{u+=Q} 



dX. 



(4.15) 



{/!+>K+}n{H+>0} L 



{h^y-{u^r 



dX < 



< 



{h+>u+}(^{u+>Q} 

{/!+>«+ }n{i(+>o} 
\h — u\'^dX. 

Br 



{h+-u+ydX 
{h-uYdX 
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Analogously, we obtain 



(4.16) 

with 
(4.17) 



Br L 



(h-y-(u-y 



dX < 



+ 



{h->ir}n{ii->0} 



{h-y-{u-y 

{u-hYdX, 



dX 



{h->ii-}r\{u->0} L 



{ir-y-iu-y 



{h->ir}n{u-=0} 



dX < I (h--u^ydX 

{h->u-}n{u->o} 



{u-hYdX 



< / \h-u\'^dX. 
Jbr 

Hence, we find 

(4.18) / Fy{h)-Fy{u)dX<C j \h-u\'^dX, 

•'Br J Br 

where C = C A_) is a positive constant. 

Combining (14.121 ). (14.1 II ) and employing Holder inequality followed by Poincare inequality 
and (14.181) we obtain 

I \V{u-h)\PdX<C3 f Fy{h) - Fy{u)dX 

JBr J Br 

<C4 [ \u-h\^dX 
Jbr 

<C5[j^ \y{u-h)\Pdxy'''\BR\'+yi"~"^ip, 



where C4 and C5 depend on p, n, and Thus, by Young inequality we reach the following 
estimate 

^/(P"^)|R„|i+iM/'r/(p-r)) 



/ Fy{h)-Fy{u)dX < C{p,y)[C{p,nM^-)y'^'' ^\b, 

J Br 



1 

4' 



(4.19) 



< c{p) [c{p,nM.^^)Y'^'''\BR\'^'i<pyi^p-y^^ 
1 

4' 



+ -jV{u-h)\\P^j„ 



where C{p, 7) = [^j '' (^V^j ^^'^ ^(p) ~ (^p ) ' • Holder inequality and Poincare inequality 
yield 

(4.20) 



4 \ 



/ f{X){h-u)dX < \\f\\L->\BR\ 
Jbr 



p-i 1 

p ''\\u — h\\LP 



< 



p-i 1 I 1 
L<i\BR\ P 9 "\\^{u-h)\\LP. 
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Thus, applying Young inequality once more, we reach 

(4.21) / f{X){h-u) < C(p)(||/||L0^|5i;|^^'^"^+"^||V(i<-/z)||£ 



I Br 

+ -\\W{u-h)\\LP 



p ] 1 1 r (g-")p ] 1 

Cip){\\fy')^\BR\ ^"^^"-'^^'U-\\V{u-h)\\LP. 



Replacing (14.191) and (14.211) in (14.101) we easily obtain 

/ \Vu-(yu)r\PdX < C(n,p,a„) - / \Vu{X)-{Vu)j,\PdX 

J Br J Br 



P 1 4- M 



+ C{n,p,ap){\\f\\L'>)p-'\BR 
< CQ"^'"' j^\Vu{X)-{VuUdX 

whereC = C(n,p, A+, A_, a^, H/Hl?) is a positive constant. In view of Lemma lZTl and W^'P bounds 
of u we conclude 

(4.22) / \Vu-{Vu)r\PdX < C(n,p,A+, A_, ||/||i,(ii),dist(Xo, ^fl)) ■ r«, 

for a entitled in (11.61) . Finally Campanato's embedding Theorem (see for instance ifTTl ) gives the 
desired Holder continuity of the gradient of u. The proof of Theorem 1 1.1 1 is complete. □ 



5 Log-Lipschitz estimates 

In this Section we address sharp regularity for jets and cavities type problems, i.e., 7=0, with 
sources in the conformal threshold case / G L"{^), where n is the dimension of the ambient. 
Hereafter u = uq denotes a minimizer of the energy functional 

(5.1) /'o(v) := (|Vvr + A+Z{v>o} +A-Z{v<o} +/(^) ■ v) dX, 

for scalars < A- < < oo. Existence and pointwise bounds for uq is has been assured by 
Theorem 13. 11 

Proof of Theorem \1.2\ We start off by fixing an arbitrary point Xq E O. and R> such that 
R < dist(Xo, do.). As before, we denote Br := Br{Xo). We follow the initial steps of the proof of 
Theorem ll.il Let h be the p-harmonic function in Br that agrees with u on the boundary, i.e., 

Aph = OmBR and /z - m G (5/;) . 



REGULARITY FOR TWO-PHASE FREE BOUNDARY PROBLEMS 20 

By Lemma l4n we have 

r / r\ n+pa,, r 

(5.2) / \Vu{X)-{Vu)fdX < C(-) / \Vu{X)-{yu)j^\PdX 

+ C [ \Vu{X)-Vh{X)\PdX. 

JBr 

On the other hand, by the minimality of u we have 

(5.3) I {\Wu\P -\Wh\P)dX < I {Fo{h)~Fo{u))dX+ f f{X){h-u)dX. 

J J B^ J B^ 

Readily one verifies that 

(5.4) / {FQ{h)-FQ{u))dX<C{X+,X-)\BR\. 

JBk 

As before, applying Holder inequality and afterwards Poincare inequality we obtain 

(5.5) / f{X){h-u)dX < \\f\\L"\BR\'l^-»\\u-h\y 
Jbr 

< \\f\\L'-\BR\'^'"+"\\V{u-h)Up. 
Therefore, with the aid of Young inequality we estimate 

/ nX){h-u)dX < C{p)i\\f\\L'')^\BR\^^'^^\\y{u-h)\\l 
Jbr 

(5.6) + h^\V{u-h)\\LP 

= C{p){\\f\\Ln)l^^\BR\+^-\\V{u-h)\\u,. 
Taking into account (15.21) and replacing (15.41) and (15.61) into (15.31) we reach 

r / r\ n+pa,, r 

/ \Vu-{Vu)r\PdX < C(n,p)(-) / \Vu{X)-{Vu)JPdX 

JBr J Br 

+ C(n,p)[C(A+,A_)]|5^|+C(n,p)C(n,p,A+,A_,||/||L«)|5^| 
/ \Vu{X)-{Vu)j^\PdX + CR'\ 

JBr 

where C = C{n,p,X+,X-, ||/||l") is a positive constant. In view of Lemma 12. 7 1 we obtain 

(5.7) / \Vu-{Vu)r\PdX <Cr^, 

JBriXo) 

which shows that the gradient of u lies in BMO space and for any fixed subdomain Q.' m Q., there 
holds 

I|V"IIbmo(£^') < c{^',n,p,x+,X-, II/IIl'O- 



r \ n+pa,, 

^ ^^r) J Br 
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From Fefferman-Stein BMO Characterization Theorem, see |[9l, there exist vector fields Fq, Fi , ■ ■ ■ F,, e 
U°{Q.'), such that 

VM(X)=Fo(X) + £^Xr;), 



i=l 



where denotes the classical Riesz transform, 



^j{f):=f*Kj for Kj{Xj):=^^. 

It now follows by a similar reasoning employed in the Appendix of lfT2l that 

\Vu{X)\<-\og\X-Xo\, forXe5p(Zo), p < 1. 
Finally, by Morrey's type estimate, we obtain, for s > n, 

\u{X)-u{Xo)\ < C\X-XQ\^--r ■ ( [ |Vm(Z)|VzV 

\JBr{Xo) J 

< C\X-Xo\^--i |logZ|'-|Z|"-ijZj 

< C\X-Xo\-\log\X-Xo\\, 

and the proof of Theorem 1 1.2l is concluded. □ 

6 Lower gradient bounds 

From this Section on, we aim towards gradient estimates to minimizers of heterogeneous p-jet 
flow functional (I5.1I ). We remark once more that even for equations with no free boundaries, say 
A = A+, it is not possible to obtain pointwise control of the gradient of uq, under the borderline 
condition f E L". In this case, as proven in Theorem II .21 the best control available is of logarithm 
order. Therefore, from this Section on, we shall assume the source function f{X), appearing in 
functional (15.11) is ^-integrable, for q> n. Under such natural hypothesis, our next Theorem shows 
that Uq grows linearly away from the free boundary : = d {u > 0} D Q.. 

Theorem 6.1. Let uq be a local minimizer to ^q, with f G L'^{Q), q> n, Q! (^Q. and Xq G {mq > 
0} n Ql. There exists a constant cq > depending only on n, p, and such that 

m(Xo) >codist(Xo,5+)- 

Proof. Let us fix Xq G {uq > 0} fl Q.'. It suffices to show such estimate for points Xq G {uq >0}nQ.' 
such that 

< dist(Xo,g'+) < A+, \\f\\ui(a)). 

for 5{n,p,X^, WlWifiia)) to be regulated a posteriori. Let us denote d := dist(Xo,^^+) and if we 
define 

v{X) :=^uo{Xo + dX), 
d 
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one easily verifies that v is a local minimizer to 

X(^):= / {\W^\P + X+X{^^^}+d-f{XQ + d-X).^{X))dX, 

mWh'^ {B\)+v. The thesis of Theorem l6.1l is equivalent to proving that v(0) is universally bounded 



away from zero. Clearly v > in 5i. By Hamack inequality (see Theorem [23]), we have 



(6.1) 



v{X) < C{n,p){v{0) + \\d-f{Xo + d-X)\\[-l^^, 



< Cin,p){viO)+{d'-'^-\\f\\,]" ■ } VXG53/5. 

In the sequel, we choose a nonnegative, smooth radially symmetric cut-off function \ff satisfying 

(^ = Oin5i/io and ^ = lin5i\5i/2 
and define the test function g mB\ by 



1 

p-i 



g{X) :=min v,C(n,p) v(0)+ J 



Notice that g E W^'^ and from Harnack inequality, estimate (16.11) . g agrees with v in Bi \ Bij2- Let 
us label the set 



Bi/2 3 n := |y G 5i/2 : C{n,p) |v(0) + [d'-'i ■ " ' | • w{Y) < v{Y) j D B^ 

From the minimality of v, we estimate 

(6.2) J^X+{l-Xis>o})+d-nXo + d-X)-[v{X)-g{X)]dX<J^{\yg\P-\yv\P)dX. 



p-1 



/lO- 



(6.3) 



n 



C{n,p){v{0) + {d'~'i 



d 'I 



p 



IV^IIc 



The right-hand side of (16.21) is readily estimated as 

{\Wg\P -\Wv\P)dX < 

< Cv{0)P + C 

We now turn our efforts towards estimating the left-hand side of (16.21) by below. Readily we obtain 

(6.4) - X{s>0}) dX = I^X+X{s=o}dX 

> A,+ \Bi/iq\. 

Livoking once more Harnack inequality (|6.1I) and the fact that n C 5i/25 we estimate 



n 



(6.5) / d-f{Xo + d-X)-[v{X)-g{X)]dX<Cld^ ? 



v{0) + {d^-'i 



1 
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Combining (lO) . (l64l) and (|631) we reach 



1 — ii 



p-1 



(6.6) C[v{0r+(d'-"^ ■ 11/11,) v(0)} > A+ |5i/io| -C 

Therefore, choosing < J < A+, ||/||L«(f2)) ^ 1^ appropriately, we conclude 

v(0)>c(n,77,A+, 11/11,) >0, 
and the proof of Theorem 16.11 follows . □ 

Next we iterate linear growth established in Theorem l6.1l as we obtain a stronger non-degeneracy 
property for uq near the free boundary. 

Theorem 6.2. Let uq be a local minimizer to ^q, with f G U^{Q), q > n, Q' Q andX^ G {mq ^ 
0} n Q!. There exists a constant > depending on n, p, and ||/||L9(fi). ^^<^h that 

sup > £ ■ r, 

S,(Xo) 

for anyO<r< dht{dQ', dO.). 

Proof. By continuity, it suffices to show uq is strongly non-degenerated, i.e., the thesis of Theorem 
16. 21 holds within the positivity set 

a+ ■= {uo>0}nQ.'. 

We will obtain such a result by iterating linear growth estimate. More precisely we will initially 
show that there exists a 5o > that depends only on n, Q.', p, A+ and ||/||, such that if X G {uq > 
0} n a', there holds 

(6.7) sup Mo > (1 + 5o)uo{Xo), 

where d{X) dist(X, ^~'"). In order to verify (16.71 ). let us assume, for the purpose of contradiction, 
that no such a do exist. If so, it would be possible to find sequences 5j = o(l) and Xj G {uq > 
0} n Q.' satisfying 

(6.8) sup Mo < (1 + 5j)uo{Xj), for dj := dist(Z;,^+) = o(l). 

BdjiXj) 

Let us consider the following normalized sequence of functions pj: Bi — )> M defined by 

_ UQ jXj + djZ) 
' •" UQiXj) ■ 

Clearly, pj{0) = 1, and from (lOl) . 

(6.9) 0< pj<\ + dj in Bi. 
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In addition, Pj satisfies 

(6.10) A,p, = -^^ipj-/(X,+J,Z), 

in the distributional sense in Bi . Taking into account the linear growth established in Theorem l6.1l 
and Equation (16.101) . we reach 

(6.11) I ApPj I < Cdj ■ fix J + djZ) , Bx. 

From Harnack inequality, we deduce the sequence {PyjyeN is locally equicontinuous in Bi \ thus, 
up to a subsequence, pj — p locally uniformly in Bi. Harnack inequality further reveals that for 
any |X| < r < 1, there holds 

(6.12) < [1 + 5j\-pj{X) < Cr {[\ + 5j\ - pj{0) +d]~'^ ■ 11/11,^ = C, -0(1). 

Letting j — oo in the above estimate, we deduce the limiting blow up function p = 1 in 5i . 

We now show that such a conclusion drives us to an inconsistency. To this end, let Yj E be 
such that dj = \Xj — Yj\. Up to subsequence, there would hold 

which clearly gives a contradiction for 7^1. We have shown the validity of estimate (16.71) . 

To finish up the proof of Theorem 16.21 we employ a Caffarelli's polygonal type of argument. 
That is, we construct a polygonal along which uq grows linearly. Starting from Xq = X, we find a 
sequence of points {X„}„>o such that: 

1. uo{Xn)>{l + 5oruo{Xo) 

2. |X„-X„_i|=dist(X„_i,5^+) 

3. uo{Xn) -uo{X„^i) > c\Xn-X„^i\. In particular, uo{Xn) - uo{Xo) > c|X„-Xo|. 

Since u{xn) °° as n oo this process must be finite, that is, there exists a last X„y in the ball 
Br{Xo). For such a last point, 

\X„^^-Xo\ > Cpr, 

Finally, 

sup Uq > Uo{Xno) > Mo(^o) + - | > C-r, 

Br{X) 

and the proof is concluded. □ 
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7 Stability for free boundary problems 

In this section we show the stability of the family of free boundary problems obtained by the 
minimization of the non-differentiable functionals 

(7.1) /y{v):= I {\Vv\P + X+{v+y + X-{v-y + f{X)-v)dX -^mm, 

J f2 

as 7= o(l). The ultimate goal of this section is to show that any limit point uq of {uy}j=Q(^i^^ is a 
minimizer to the /^-degenerate cavitation functional 

(7.2) ^o(v) := (|Vvr + A+Z{v>o} +^-Z{v<o} +/(^) ■ v) dX. 

Initially we show compactness of {My}o<y<i in the W^'P topology. 
Proposition 7.1. Let Ujj be a sequence of minima to the functional ^y., f G L" and assume Uy. — v 

1 E 

a.e., Yj 7o- Then for any < E < oo, Uyj ^ v in the WjJ^ {Q) topology. 

Proof. It follows from Proposition 13.31 and a.e. convergence that ApUyj ApV in the sense of 
measures. Thus, from truncation arguments, see for instance BH, 

(7.3) Vuyj Vv a.e. in D.. 

From Theorem 1 1.21 for any Q.' d Q., there exists a constant C(n,p, A+, A-,n', H/H/O^ independent 
of Yj, such that. 



(7.4) llVMyjBMO(n') <C{n,p,X+,X-,Q.', 
Thus, from John-Nirenberg's Theorem, for 1 < £ < oo fixed, 

(7.5) ||VMy.||^£+i(fj,)C(n,/?,A+,A_,f2', ||/||„). 

Finally combining (17.31) . (17.51) and classical arguments, see for instance, lfT3l . we deduce 

Vwy. ^Vv in L^{Q.'), 
and the Proposition follows. □ 
Theorem 7.2. Let uq := limwy^. as Yj 0. Then uq is a local minimizer of ^o- 

Proof. Let B, be a ball in Q.. Given an arbitrary W^'^ function y that agrees with uq on dBr, we 
have to show that 

By density we may further assume that i// is bounded. Let us define the interpolated function 

{\X\ — r 
uo + —^{uyj-uo) in B,+h\B, 
1// in Br- 
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One simply verifies that 



(7.6) \yxiry^,h\P < Cp!^\Vuo\P + ^\uy^-uo\P + \yuy^-Vuo\pY mB,+h\B,. 



In the above estimate, we have used the classical facts: 



(7.7) 



V(|Z|) = ^ and 



\X\-r 
h 



<1 in B,+h\Br. 



By L°° bounds, Theorem 13. 1[ there exists a constant Ci > 0, independent of Yj, such that ||My^.||oo < 
Ci. Thus, if we denote 

Hp.it) := (r±)r., 

we have 
(7.8) 
and 

(V.9) Hfji¥rj,^<{\ML-iB,)y'X{u,.m^ 
We can estimate 

^yj{B,+h,\iry.M) = [ \y\iryj,h\P + ?i+H+{\iry.j,)+?i^H^.{\iry.j,)dX 



(7.10) 



+ / ^ f{X) 



< 



UQ H [Uyj - Uq) 



dX + /j.{Br,W) 



Cp I \Vuo\ PdX +Cp I I Vuyj - Vmo I PdX 

+ [2X+{3Ciyi +3Cx\Cp\Br+h\Br\ + % f \Uy.-UQ\PdX 

+ ^o(5.,V/) + (||v/||;i(g^.)-l) / A+Z{v.>o}+^-Z{v/<o}^^ 

+ \Br+h\Br\ '>\\f\\ui{Q.)- 

By pointwise convergence Uyj uq we have 

(7.11) lim / luy. — uol^dx — 

and by Proposition l7.1l 



(7.12) 



lim / \Vuy - Vuo\Pdx = 0. 

j^'-'JBr+h\Br 
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From the minimality property of Ujj, 

(7.13) fYj{Br+h. ¥rj,h) > /Yj{Br+h.Uyj) > ^y.{Br,Uyj). 

Furthermore, it follows from Proposition |7T| 

(7.14) / \Vuo\PdX=\im [ \Vuy.\PdX. 

J Br j^'=°JBr 

By the pointwise convergence Uy^ — )■ uq and Fatou's Lemma (see the proof of Theorem 13.11) . we 
conclude 

(V.15) ^ KX{ua>Q} + ^-X{u,,<Q}dX < limi;^f^ K ("r,)^'^^r,>o} {^yi)^^ X{uy.<Q}dX, 
and 

(7.16) lim / f{X)uy.dX=\\m [ f{X)uodX. 

Finally, combining (I7.10I) - (I7.16I) we reach 

^o(5,-,wo) < liminf^y (5^+;„My ) 

J— >«j 

(7.17) < MBn¥)+Cp[ \Vuo\PdX 

+ (2A+ + 3Ci)Cp\B,+h \ B,.\ + \B,+h \ B,\ . 
Letting /j — 0, we finish the proof of Theorem r7!2l □ 

8 Existence of Lipschitz minimizers for 

In this final Section we conclude the proof of Theorem II. 3 [ 

Proof of Theorem [O] We will show that any limit point of a family of local minimizers of ^y, 
as 7 — ;> is locally uniformly Lipschitz continuous. Theorem 1 1.3 1 thus follows from this fact and 
the stability property proven in Theorem l7!2l 

Initially, we point out that from the estimate obtained in Theorem 11.11 and classical Cam- 
panato's inequality, it suffices to show the existence of a interior point E ^ such that 

(8.1) \yuy{Z,)\ < C(n,p,A+, A_, 11/11,) < +00, 

for a constant C(n,p,A+, A- , \\f\\q) that is independent of 7^ 1. 

Our first observation is that, since with no loss of generality, we are under the situation where 
Uq ^ 0, there must exist a point E Q. and a small constant 5-k > 0, such that 

(8.2) dist(Z^,(9{My>0}) > 5,, >0, V7G [0,1]. 
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By nondegeneracy of uq, Theorem l6.1[ there holds 

(8.3) uo{Z^)>c-5^. 

For < r < 5^, we apply Hamack inequality, Theorem l2.3[ to conclude that 

(8.4) sjfl.)"" - '^^^{^■^^'i''^' ■WfW^w)'"'] 

> c^, 

provided 

(8.5) < r < r A+, ||/||^) =: < 1. 
Since Uy converges locally uniformly to uq, we have 

(8.6) inf Uy > ^Ci,, 

provided 7< 7* ^ 1. From Proposition 13. 3[ Uy solves, in Br^{Zi,), the following equation, in the 
distributional sense, 

(8.7) ApUy=h+{u+y-' + ^nX) =: g{X). 

Recall / G L^, for q> n, thus, from (18.61) . g EL^ {Br^{Zi,)), for ^ > n. By Local Lipschitz regularity 
for heterogeneous equations, see Theorem 1.1 from Q, we finally estimate 

(8.8) |VMy(Z^)| <C{n,p,\\uy\\oo,r^,c^,\\f\\q) 

and the proof of Theorem 11.31 follows . □ 
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